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5_i ■ ABSTRACT 

We present a fully sampled C ls O (1-0) map towards the southern giant molecular cloud 
(GMC) associated with the Hll region RCW 106, and use it in combination with previous 
13 CO (1-0) mapping to estimate the gas column density as a function of position and velocity. 
We find localized regions of significant 13 CO optical depth in the northern part of the cloud, 
OO ' with several of the high-opacity clouds in this region likely associated with a limb-brightened 

0^ \ shell around the H II region G333.6— 0.2. Optical depth corrections broaden the distribution of 

CnJ ■ column densities in the cloud, yielding a log-normal distribution as predicted by simulations of 

CS| ' turbulence. Decomposing the 13 CO and C ls O data cubes into clumps, we find relatively weak 

, correlations between size and linewidth, and a more sensitive dependence of luminosity on 

■ size than would be predicted by a constant average column density. The clump mass spectrum 
CyQ \ has a slope near —1.7, consistent with previous studies. The most massive clumps appear to 

have gravitational binding energies well in excess of virial equilibrium; we discuss possible 
explanations, which include magnetic support and neglect of time-varying surface terms in the 
virial theorem. Unlike molecular clouds as a whole, the clumps within the RCW 106 GMC, 

■ while elongated, appear to show random orientations with respect to the Galactic plane. 

;_j 

C$ ' Key words: ISM: molecules - ISM: clouds - ISM: structure - stars: formation 



1 INTRODUCTION 

Molecular clouds must be supported by turbulent motions in 
order to resist gravitational collapse and thereby explain the 
relatively long time-scale for star formation in the Galaxy 
dZuckerman & Palmejl974) . Observations of clumpy and filamen- 
tary structures in these clouds, in qualitative agreement with numer- 
ical simulations of turbulence, provide further evidence for turbu- 
lent support. However, the simulations also indicate tha t turbulence 
dissip ates quickly (within a few crossing times, e.g. IStone et all 
1998), and must therefore be continually renewed by energy in- 
jection. Given that several physical processes, acting on different 
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scales and contributing different amounts of kinetic energy, can 
contribute to turbulence, the nature of the principal driving source 
of turbulence in molecular clouds remains a subject of active de- 
bate. 

When studying the turbulent properties of molecular clouds 
and their possible universality, statistics which distinguish among 
different spatial scales are generally pre ferred: these include 
Fouri er power spectra of 2-D velocity slices dLazarian & Pogosvanl 
200G), the delta variance ( Stutzki et al. 1998), and principal compo- 
nent analysis faever & Schloerbll 19971) . These techniques take ad- 
vantage of the self-similar nature of turbulence to reduce its charac- 
terization to a small number of quantities, e.g. power-law indices of 
the energy spectrum or of the correlation between spatial scale and 
velocity dispersion. A number of studies in recent years have ap- 
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plied these methods to atomic and molecular line data I Dickey et al.l 
1200 ll ; iBensch et alfeOOll ; iBrunt & Heverll2002l ; ISun et alj2006h . 

On the other hand, relating turbulence to star formation re- 
quires that the molecular gas data be considered as a function 
of position, and not only as a function of spatial scale. To facil- 
itate such analysis, it is commonplace to decompose molecular 
line maps into discrete structures, usually referred to as clumps or 
cores, and to relate the properties of these clumps and cores to the 
properties of stars and star-forming units. Although several auto- 
mated deco mposition methods have been introduced (e.g., GAUSS- 
CLUM PS, lstutzki & Guestenll990l : CLUMPFIND. IWilliams et al.l 
all are sensitive to the map resolution and to ad hoc param- 
eters selected by the user. Papers emphasizing numerical results 
have t ended to cast doubt on the interpretation of these decom posi- 
tions ( iBallesteros-Paredes & Mac Lowll2002l ; iDib et alj|2007h . yet 
they remain useful as a way to identify regions that are amenable 
to multi-wavelength analysis. 

Motivated by an interest in both approaches to studying 
molecular cloud structure, we have recently undertaken a system- 
atic program of observing the giant molecular cloud (GMC) asso- 
ciated with the RCW 106 star-formin g region in multip le molec- 
ular lines with the Mopra Telescope. iBains et alj hood hereafter 
Paper I) presented the initial results of the program. Applying a 
two-dimensional dump-finding algorithm to an integrated 13 CO 
( J=l — *• 0) emission map, we found that about half of the identified 
13 CO clumps could not be d istinguished as clumps in the 1.2-mm 
dust emission map made by Mookerjea et al.1 J2004). While some 
of these clumps may have fallen below the sensitivity limit of the 
dust map, it is possible that a significant fraction of the gas clumps 
were the result of integrating unrelated cloud structures in velocity. 

In this paper we present our C 18 ( J=l — > 0) mapping of this 
region, and derive improved estimates of molecular column density 
across the cloud. We also quantify and discuss variations in 13 CO 
optical depth and their implications for estimating the probability 
density function (PDF) of the column density. Finally, we revisit the 
structural decomposition of the cloud emission, working directly 
on the three- dimensional data cubes, using the CPROPS algorithm 
developed bv lRosolowskv & Lerovl d2006 , hereafter RL06). Future 
papers in this series will examine the power spectra of the molecu- 
lar gas (P. Jones et al., in preparation) and the relationship between 
recent star formation and the molecular clumps (N. Lo et al., in 
prep aration). As in P aper I, we adopt a distance to the cloud of 3.6 
kpc jLockmanl 1 1 9791) . giving it a projected size of roughly 90 x 
30 pc. Since all observed lines are J=l — * transitions, we simply 
refer to them hereafter by the name of the species (' 13 CO', etc.). 



2 OBSERVATIONS 

Observations of the C ls O (J=l — > 0) line towards RCW 106 were 
made with the 22-m Mopra telescope of the Australia Telescope 
National Facility (ATNF) between 2005 July 5 and September 1 1 . 
A total of 46 5' x 5' fields were observed in two passes, with the 
first pass produced by scanning along the R. A. axis and the second 
pass along the Dec. axis. The fields observed in each pass were cho- 
sen to have some overlap, with centre positions separated by A'.15, 
as for the 13 CO observations described in Paper I. In addition, the 
field centres for the second pass were shifted by 1' in both R.A. 
and Dec. (towards the southwest) in order to further reduce sensi- 
tivity variations at the field edges. The region mapped is shown in 
Figure[T] and is a subset of the region previously mapped in 13 CO, 
which comprised 93 fields. 
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Figure 1. C 18 map integrated from —65 to —35 km s -1 . Contour levels 
are linearly spaced starting at 7 K km s^ 1 (T m b). The dashed line repre- 
sents the region mapped with Mopra, while the larger solid line represents 
the region previously mapped in 13 CO (Paper I). The Mopra beam size is 
shown in the lower left corner. The orientation of the Galactic coordinate 
axes is also shown at lower left. 



The observations were conducted in raster scan ('on-the-fly') 
mode, with the same mapping parameters and speed as for 13 CO 
(Paper I). Each field took ~70 min to complete, and the system 
temperature was calibrated every ~30 min using an ambient load. 
Typical above-atmosphere system temperatures were ~300 K. The 
pointing solution was refined after each field using observations of 
the nearby SiO maser IRSV 1540; typical corrections were ~5". 

The Mopra digital correlator was configured to output 1024 
channels across a 64 MHz bandwidth, providing a velocity reso- 
lution of 0.17 km s _1 . To ensure that reference and source obser- 
vations were taken at the same frequency, Doppler tracking was 
disabled and the center frequency was set to 109.792 GHz, which 
corresponded to vlsr ~ —50 km s" 1 at the time of the observa- 
tions. As in Paper I, we quote velocities in this paper relative to the 
local standard of rest (LSR). 

Data reduction was performed using the ATNF MIRIAD, 
LIVEDATA and GRIDZILLA packages^ Following data in- 
spection and editing, spectra were bandpass corrected using an 
emission-free spectrum taken at the beginning of each row at a fixed 
sky position (a2ooo=16:27, ^2000=— 51:30). A linear baseline was 
fitted to channels deemed to be emission-free and subtracted. The 
spectra were then gridded into a datacube using a 33" (FWHM) 
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Figure 2. Average C ls O spectrum, in T mb units, over the region shown in 
Fig.[T] The thick line is the 13 CO spectrum averaged over the same region, 
divided by the adopted abundance ratio of 7.4. 



Gaussian smoothing kernel and 15" pixel scale. Approximating the 
Mopra point source response as a Gaussian with FWHM 33", the 
final map resolution was determined to be 45". The RMS noise 
in the map was 0.12±0.07 K (T£) and ranged from 0.063 to 1 K; 
high-noise regions at the edge of the map were masked out by set- 
ting a threshold RMS of 0.25 K. The resulting masked map has an 
RMS noise of 0. 1 1±0.03 K. FigureQ]shows the C 18 intensity in- 
tegrated from —65 to —35 km s _1 , a velocity range which covers 
the bulk of the emission. 

Since the Mopra beam response is not perfectly Gaussian, a 
correction factor is needed to account f or flux which is s cattered out 
of the main lobe of the primary beam. lLaddet all (2005) derived a 
main-beam efficiency for the 2004 season of rj m b=0.42 at 1 15 GHz, 
and an 'extended' beam efficiency (for sources that couple to the 
inner error beam) of 77^=0.55. Because the structures of interest 
for this paper are comparable to or larger than the size of the inner 
error beam (~80"), the extended beam efficiency is used here. To 
determine whether there was a change in r] x b from 2004 to 2005, we 
compared uncorrected (TJ scale) spectra taken by one of us (BM) 
in 2004 June towards 15 positions in the RCW 106 cloud with the 
2005 map data (with minimal smoothing applied in the gridding 
process). For 449 channels with T£>1 K, the average ratio of 2005 
to 2004 intensities was 0.9 with a standard deviation of 0.18. Since 
the data suggest a change of ~10 per cent in the efficiency, we 
adopt a value of 7^6=0.5 to correct the intensity scale to T m b units, 
although the intensity calibration is likely to remain uncertain at the 
~ 10-20 per cent level. 

Figure|2]shows a spectrum of the C 18 brightness temperature 
averaged over the map. The 13 CO brightness temperature over the 
same region is also shown, scaled down by the assumed abundance 
ratio of 7.4 (see ij3.lt . Aside from calibration errors, any mismatch 
between the profiles could be due to variations in the abundance 
ratio and/or to optical depth effects. 



3 ANALYSIS 

3.1 Opacity, excitation temperature and column density 

Column densities were derived from the 13 CO and C ls O data by 
assuming a filled beam and a uniform excitation temperature, com- 
mon to both tracers, within the beam. The first step was to derive 



the opacity in both lines at each (a, 5, v) volume pixel (voxel) in 
the datacube from the ratio of 13 CO to C ls O main-beam brightness 
temperatures (e.g. jMvers et aljl983h : 



r m b( 13 CO) _ 1 - exp(-ns) _ 1 - exp(-7.4r 18 ) 



T mb (C 18 0) 1 - exp(- 



•TlsJ 



1 - exp(-Tig) 



(1) 



To reduce the effects of noise, the spectra were binned to a velocity 
resolution of 0.68 km s _1 , and the brightness ratio was only com- 
puted when the signal-to-noise ratio (SNR) was >5 in the C 18 
data. Voxels which failed this test, or had brightness ratios <1 (all 
such voxels were found to be consistent with noise), were assigned 
13 CO column densities in the optically thin limit. 

Equation Q] assumes ri3 = 7.4ri8, based on the abundance 
ratio [ 13 Clf 16 Ql/[ 12 Ci r 18 Ql=7.4 estimated a t a galactocentric ra- 
dius of D G c=5.5 kpc Jwilson&Roodlll994l) . Note that this value 
is ~30 per cent highe r than the value of [ 13 CO]/[C ls O]=5.5 used 
bv lMvers etai] fl983), which assumes terrestrial abundances. Our 
estimate, while still quite uncertain, takes into consideration the 
isotopic abundance gradients in the Galaxy and the decrease in 
[ 12 C]/[ 13 C] since the epoch of the Earth's formation due to CNO 
processing. 

The excitation temperature T GX can then be derived from the 
equation of radiative transfer. 



Tmb = /[J(Tex) - J(T bg )][l - exp(-r)] 



(2) 



where / is the beam filling factor, J(T) = T /[exp(7o/T) - 1], 
T bg =2J K, and T = hv/k = 5.29 K for the 1-0 transition of 
13 CO. This equation was applied to the 13 CO rather than C ls O 
data because of the higher signal-to-noise of the former. For most 
of the analysis we assume a filled beam, so f=l. However, in cases 
where the emission was weak but r was found to be significant, re- 
sulting in T cx <5 K, we set T ox =5 K and calculated the filling factor 
/ needed to account for the observed brightness temperature. Such 
regions of very low implied T cx tend to be found at the edges of 
the cloud (in position and velocity), where incomplete beam fill- 
ing is quite likely. We also imposed a maximum T13 of 12, corre- 
sponding to a brightness temperature ratio of 1.25, as we cannot ac- 
curately measure higher opacities given uncertainties in telescope 
calibration and the abundance ratio. Such voxels were extremely 
rare (<0.02 per cent) and only found in the highest opacity region 
in the northwest part of the map. 

Given T cx , T13, and Tig, the column densities of 13 CO and 
C ls O can be derived using the expression for column density in a 
single rotational level (e.g. jGarden et all 199 lh : 



Nj = 



3h 2 J + 1 



8tt 3 /j 2 J + 1 



exp 



(-—) 



dv 



(3) 



where J is the lower state of the transition and fi is the perma- 
nent dipole moment of the molecule (n=0. 1 1 for the lines studied 
here). A lower limit to the total column density Ntot can be derived 
by assuming only th e first two levels (J=0 and 1) are populated 
I IFrerking et alJI 19821) . while an upper limit can be derived assum- 
ing all rotational levels are thermalized with the same excitatio n 
temperature T cx . The latter assumption yields jGarden et al.ll99lh : 



3A- 



T cx + hB/3k 



dv 



(4) 



8tt 3 B/j 2 1 - exp(-hu/kT eK 

where B is the rotational constant. 

In a significant fraction of voxels (roughly 12 per cent of the 
total), generally with detected but weak C 18 emission, we found 
T mb ( 13 CO)/T mb (C 18 0) > 7.4, resulting in optical depths which 
were undefined in EquationQ] The largest significant values for the 
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Figure 3. Images of integrated column density and peak opacity. Left: Integrated H2 column density assuming all rotational levels are thermalized. Contour 
levels are 6.4n 2 X 10 21 cm 2 (n=l,2,3,4,5). This gives an upper limit to the column density. Centre: Integrated H2 column density assuming only the J=0 
and 1 states are populated. Contour levels are again 6.4n 2 X 10 21 cm 2 . This represents a lower limit to the column density. Right: Peak 13 CO optical depth 
at each position in the map. Contour levels are 2 An (n=l, 2,3,4,5). 



© 0000 RAS, MNRAS 000, 000-000 



Mapping the column density in the GMC associated with RCW 106 5 



ratio are ~14 and occur near (a, S, v) = (16:20:53, —50:43:55, 
—56). These high values could result from a slight mismatch in 
beam sizes near regions of bright emission, from differences in 
excitation between 13 CO and C ls O, or from chemical fractiona- 



ising the al 

dLanger et al.ll98fj|) . For these cases we assumed an excitation tem- 
perature based on the 12 CO brightness temperature as described in 
$3.2[ and derived iV(C 18 0) in the optically thin limit, where 



T ex + 0.88 



T mb » [J(T«) - J(T bg )]r . 
This yields (cf. lBourke etai1ll997h : 

7V(C 18 0)thin = 2.42 x 10 1 



1 - exp(-5.27/T ox ) 

1 / - ,„18 



J(T CX ) - J(T bg ) 



T mb (C ls O) dv . (5) 



We also applied this procedure whenever Ti3<0.1, to limit the 
number of poorly constrained, very high values of T cx . If the col- 
umn densities had been calculated from 13 CO instead, they would 
have been on average ~10 per cent higher (assuming an abundance 
ratio of 7.4), though in some regions up to a factor of 2 higher. 

In other voxels, only 13 CO was observed or detected. Again, 
an excitation temperature was assumed based o n 12 CO, and 
N( 13 CO) was derived in the optically thin limit dBourke et al.l 
E993): 



7V( 13 CO)thm = 2.42 x 10 1 



T cx + o.; 



exp(-5.29/r ex ) 



J(T ex ) - J(T( 



hg ) 



T mb ( 13 CO) dv . (6) 



The molecular hydrogen column density, N(H.2), was then 
calculated by integrating the 13 CO or C 18 column density cube in 
velocity from —80 and — 20 km s" 1 assu ming an [H2/ 13 CO] abun- 
dance ratio of 7 x 10 s dFrerking et all 1 19821) and an [H 2 /C 18 0] 
abundance ratio of 5 x 10 6 (for consistency with our adopted 
[ 13 CO/C ls O] abundance ratio). Figure [3] shows the resulting col- 
umn density maps (assuming thermalisation of all J levels or only 
population of J=0 and 1) as well as a map showing the peak op- 
tical depth of 13 CO at each spatial position. Note that the column 
density outside the region mapped in C 18 is more uncertain due 
to our lack of constraints on the opacity. 



3.2 Default T cx using NANTEN data 

In order to provide an independent estimate of the excitation tem- 
perature in cases where it could not be derived from the 13 CO and 
C ls O data directly, we made use of the 12 CO data from the NAN- 
TEN Galactic Plane Surve y. The observat i ons an d data reduction 
have been described by iMatsunaga et alj |2001). In brief, spec- 
tra taken with a grid spacing of 4' using the 4-m NANTEN tele- 
scope (half-power beamwidth 2'6). An acousto-optical spectrome- 
ter (AOS) provided velocity coverage of 650 km s _1 at a resolution 
of 0.65 km s _1 . The resulting map was interpolated onto the Mo- 
pra grid and Tex for each voxel was derived from the peak T mb 
at the corresponding position in the 12 CO data cube, assuming a 
filled beam and optically thick CO (ri2 3> 1). Where the peak 
T mb ( 12 CO) was less than 5 K, the minimum T cx =5 K was adopted 
as above. No attempt was made to deconvolve the 12 CO map to ap- 
proximate the resolution of the Mopra data, so this method is likely 
to yield an underestimate of T cx in regions of bright emission while 
an overestimate in adjacent regions. In addition, the assumption of 



optically thick 12 CO in a filled beam is likely to yield too low an 
estimate of T cx , especially at the edges of the cloud, although our 
imposition of a minimum T cx should help to reduce errors intro- 
duced by this assumption. 

3.3 Decomposition into clumps 

Paper I had noted several velocity components in the 13 CO cube, 
with most of the emission centered on an LSR velocity near —50 
kms -1 (Fig. [2}. In order to distinguish the main GMC emission 
from unrelated emission along the line of sight, we applied the 'di- 
lated mask' method in the IDL-based CPROPS package of RL06. 
The 13 CO spectra were binned to 0.34 km s _1 resolution to in- 
crease the signal-to-noise ratio, and a mask was created by starting 
from a threshold contour of T% — 5 K and expanding in both po- 
sition and velocity to include all simply connected regions above 
a brightness threshold of 0.26 K (3<r). The resulting "GMC mask" 
was applied to the 13 CO cube, yielding a single continuous region 
encompassing 92 per cent of the emission between —80 and —20 
kms -1 and extending in velocity from —68 to —25 kms -1 . Of 
course, such a procedure does not guarantee that unrelated emission 
is not projected into the same velocity range. However, given that 
the region of interest is located somewhat away from the Galactic 
plane (from 6=— 0?2 to — 0?6), it seems reasonable to assume that 
the emission within the mask comes from a single cloud. 

To decompose the GMC into discrete structures, we applied 
a modified version of the CPROPS algorithm to the masked data 
cube. Significant peaks (local maxima) in the emission were iden- 
tified as voxels which are greater than all their neighbours within 
an I x I x k box, where 1=5' and k=3 km s _1 . The adopted val- 
ues of I and k were chosen to reduce blending between clumps and 
ensure that the decomposition could be run in a reasonable amount 
of time. However, they may bias our analysis against detection of 
closely spaced clumps. The data were then contoured below each 
peak until the highest contour level which begins to enclose an- 
other peak; this level (which is unique to each peak) is termed the 
merge level. Peaks which did not lie at least 2<r above the merge 
level, or did not contain at least 2 beam areas of pixels above the 
merge level, were discarded. The resulting 'clumps' were defined 
as the regions uniquely associated with each peak (i.e., above the 
merge level). Unlike the original CPROPS procedure of RL06, we 
did not filter the list of clumps based on the change in moments 
when merging with adjacent clumps. Such filtering, which has the 
effect of merging clumps which blend smoothly into each other, 
drastically reduces the number of clumps (from ~600 to ~30), and 
in our case leaves a great deal of substructure within each clump. 

The clumps derived from CPROPS contain a relatively small 
fraction (~10 per cent) of the total 13 CO emission, since con- 
tour levels which enclose several peaks (blended contours) are 
not decomposed. As an alternative, we also performed decompo- 
sitions using the /ECLUM P option in CPROPS, which applies the 
CLUMPFIND algorithm of Willi ams etall Jl994h to the previously 
identified emission peaks. In this case, blended contours, rather 
than being discarded, are assigned to existing clumps via a 'friends- 
of-friends' algorithm. This has the effect of assigning all of the flux 
within the GMC mask into clumps, and thus producing a 'patch- 
work' of much larger clumps. Unlike CLUMPFIND, however, the 
number of clumps cannot increase as one works down to lower con- 
tour levels: the initial number of peaks remains the final number of 
clumps. 

We also ran both procedures discussed above on two addi- 
tional cubes: a non-GMC 13 CO cube, which consists of the emis- 
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Figure 4. Locations and sizes of the clumps derived from the CPROPS analysis. Sensitivity and resolution corrections have not been applied. The ellipse major 
axis represents the clump major axis as derived from principal component analysis, (a) Integrated 13 CO emission from the GMC. Contour levels are spaced by 
10 K km s (T^). (b) Non-GMC integrated 13 CO emission. Contour levels are spaced by 3 K km s (T^). (c) Integrated C 18 emission. Contour levels 
are spaced by 3 K km s _1 (T£). 
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Table 1. Estimates of cloud mass. 





13 CO map region 


C 18 map region 


Optically thin 13 CO, all J 


2.0 


1.3 


Optically thin C ls O, all J 




1.5 


Opacity corrected 13 CO, J=0,1 


0.94 


0.69 


Opacity corrected 13 C0, all J 


2.2 


1.6 



All masses are given in units of 10 6 Mg. 



sion between —80 and —20 km s -1 that lies outside the GMC 
mask, and the whole C ls O cube within the same velocity range. 
In both cubes, islands of discrete emission were identified by start- 
ing from a threshold of 4a in two adjacent channels and expanding 
to a 3cr boundary (still requiring significant emission in at least 2 
adjacent channels). Each island was then decomposed into clumps 
as was done for the GMC 13 CO emission. Although the C ls O cube 
was decomposed without the use of the GMC mask, it turned out 
that all but two of the resulting clumps were within the mask bound- 
ary, so effectively all of the decomposed C ls O emission is associ- 
ated with the GMC. The non-GMC 13 CO cube, on the other hand, 
contains emission which appears separated from the GMC in posi- 
tion or velocity, and thus may be foreground or background emis- 
sion. 

The final step in the CPROPS analysis is to determine basic 
properties (radius, linewidth, and luminosity) for each clump. The 
radius of a clump R is proportional to oy, the geometric mean of 
the RMS sizes (second spatial moments) of the clump along the 
major and minor axes, with the major axis determined via princi- 
pal component analysi s. Following RL06, w e take R = 1.91oy (as 
previously adopted bv lSolomon et al.l 19871) . The FWHM linewidth 
AV is equal to 2.35c,,, where a v is the intensity-weighted second 
moment along the velocity axis. The luminosity is based on the 
integrated flux of the cloud and the adopted distance (3.6 kpc, ex- 
cept for the non-GMC clumps where the near kinematic distance 
was used). RL06 detail procedures to correct for the sensitivity and 
resolution biases that affect real observational data; this involves 
extrapolating the major and minor axes, linewidth, and flux to infi- 
nite sensitivity based on their values at various signal-to-noise lev- 
els, and then (for the radius and linewidth) subtracting the beam 
or channel width in quadrature. In the end we did not apply these 
procedures, for reasons discussed further in ^4. 3 1 The locations and 
sizes of the clumps are shown overlaid on the respective integrated 
intensity maps in Figure|4] 



4 RESULTS 

4.1 Overall cloud properties 

The total gas mass can be estimated from the 13 CO and C ls O cubes 
separately by assuming optically thin emission and a single exci- 
tation temperature for all rotational levels. We integrated all ob- 
served emission from —80 to —20 km s _1 using the abundance 
ratios given at the end of i]3.1l Adopting T cx =20 K as in Paper I, 
and including a factor of 1.36 for helium, the resulting masses are 
shown in the first two lines of TableQ] Note that the 13 CO estimate 
is corrected to be a factor of ~2 higher than in Paper I, due to an 
error in the original calculation. Although use of 13 CO alone may 
overestimate the column density in some optically thin regions, per- 
haps due to fractionation ( i]3.1t . it underestimates the column den- 
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Figure 5. Histogram of normalised column densities in the RCW 106 GMC, 
integrated from —80 to —20 km s _1 over the region observed in 13 CO. (a) 
Column density based on the optically thin assumption for 13 CO and adopt- 
ing T cx =20 K. (b) Column density calculated based on the optical depth 
analysis. In regions where only 13 CO is observed, N is calculated in the 
optically thin limit using T ox based on the NANTEN data. For both pan- 
els, a best-fit Gaussian has been overplotted which represents a log-normal 
distribution with the dispersion indicated. 



sity in optically thick regions, resulting in an overall mass estimate 
lower than that from C ls O. 

The last two lines of TableQ]lists the opacity-corrected masses 
derived by integration of the column density cube. The opacity- 
corrected masses, assuming population of all J levels, are only 
slightly larger than the masses derived in the optically thin approxi- 
mation, which may reflect both optical depth effects and variations 
in T ex that reduce the gas mass in cold regions and increase it in 
warm regions. 

The distribution of integrated column density values over the 
region mapped in CO, after binning to 45" pixels to approximate 
the spatial resolution of the data, is shown in Figure [5] Column 
densities < 2 x 10 21 cm~ 2 have been excluded, reflecting a 4<r 
sensitivity limit of 3 K km s" 1 (T m b) in the 13 CO data. However, 
imposing this cutoff has little effect, as <0.5 per cent of the pix- 
els fall belo w this limit. To facilitate compariso n with theoretical 
models (e.g.. lVazquez-Semadeni & Garcfal200lh . we have normal- 
ized TV by its mean value and plotted it on a logarithmic scale. The 
optically-thin, 13 CO-based column densities (upper panel) show a 
distinctly asymmetric distribution with a steep decline at high col- 
umn densities. The opacity-corrected values (lower panel) show a 
less abrupt decline towards high column densities, and the overall 
distribution can be well described as a log-normal function with a 
dispersion of 0.33 dex. Underscoring this point, the best-fit Gaus- 
sian shows a significant improvement in reduced \ 2 (assuming 
each bin has an uncertainty of yN), from 12.7 to 3.3. Thus, al- 
though the opacity correction has little effect on the total mass es- 
timate, it has a pronounced effect on the distribution of column 
densities, bringing them closer to a log-normal distribution. 
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Figure 6. Locations, sizes, and orientations of 13 CO clumps that overlap 
with C 18 clumps in both position and velocity (open ellipses) and C 18 
clumps that overlap with 13 CO clumps (shaded ellipses). The outer contour 
indicates the region mapped in C 18 0. 



Table 2. Number of clumps found with CPROPS. 





13 CO, GMC 


13 CO, non-GMC 


C 18 0, GMC 


Initial number 


594 


121 


151 


After edge culling 


456 (345) 


83 (71) 


138 (116) 


After all culling 


455 (345) 


83 (71) 


135 (114) 



Numbers in parentheses correspond to the ECLUMP method. 



4.2 Clump orientations and masses 

As noted in i )3.3l we decomposed the cloud into clumps using two 
methods: the 'default' method in CPROPS and the 'ECLUMP' 
method which applies a modified CLUMPFIND algorithm. Al- 
though the initial number of clumps from the two methods was 
the same, we then proceeded to exclude clumps which extended to 
the edge of the observed field, as their sizes and fluxes could not 
be reliably determined. Because the ECLUMP method produces 
larger clumps, it yielded more clumps which had to be excluded. 
Additional culling was then applied to remove clumps with unde- 
fined linewidths or radii, or which lay outside the GMC boundary 
in the case of the C 18 cube. As summarized in Table [2] the final 
numbers of clumps were 455, 83, and 135 in the GMC-only 13 CO 




Axis ratio 

Figure 7. Distribution of position angles (measured from north towards east 
on the sky) and minor-to-major axis ratios (b/a) for the 13 CO clumps. The 
Galactic plane has a position angle of approximately 45° (vertical dashed 
line). 

cube, non- GMC 13 CO cube, and C 18 cube respectively. We find 
a large fraction (~75%) of the C ls O clumps overlap in position- 
velocity space with 13 CO clumps, although the reverse is not true, 
as a result of the smaller region observed in C 18 and the higher 
signal-to-noise of the 13 CO observations. The locations, sizes, and 
orientations of the 223 overlapping 13 CO and C 18 clumps are 
shown in Figure [6] indicating that when there is overlap one typi- 
cally finds good agreement in terms of clump size and orientation. 

Figure |7] shows the distributions of position angles and axis 
ratios for the GMC-only 13 CO clumps, as derived using the default 
CPROPS method. In general the clumps show significant elliptic- 
ity , with a mean axis ratio ~0.6, similar to the mean value found 
bv lMvers et alj JT99ll) for cores within dark clouds. Since the axis 
lengths have not been corrected for beam smearing, the true elliptic- 
ity could be even greater. Given the filamentary appearance of the 
GMC, this is not unexpected. However the position angles of the 
clumps appear to be randomly oriented, with no clear preference 
to be aligned with the Galactic plane. We discuss this observation 
further in $5A\ 

CPROPS provides two estimates of clump mass, a 'luminous 
mass' Mlum which is simply proportional to the observed line 
flux, and a 'virial mass' Mvir defined as the mass of a spherical 
cloud with a density profile p oc r~ in equilibrium betw een inter- 
nal velocity dispersion and gravity dSolomon et aill 19871) : 

Mvir = 189 M AV 2 R 
= 1050 M @ olR , 

where a v is the velocity dispersion in km s" 1 and R is the effective 
radius in pc (R — 1.9 lay). We also calculated a 'local thermody- 
namic equilibrium mass' Mlte for each clump, being the mass in 
the corresponding voxels of the H2 column density cube assum- 
ing a thermalised population of J-levels. As Figure[8]indicates, the 
relation between Mlte and Mlum is roughly linear, with only 
a slight non-linearity appearing at larger masses, where the lumi- 
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Figure 9. Distribution of opacity-corrected masses for the clumps. The 
shaded histogram shows results from the default CPROPS analysis, while 
the dot-dashed line shows results from the ECLUMP method. The vertical 
dashed line is our adopted completeness limit. A fit to the shaded histogram 
is shown, with the corresponding value of j3. 



nous mass (essentially an optically thin estimate) underestimates 
the LTE mass by a factor of ~2.5. The virial mass, on the other 
hand, is likely to be a poor estimator of the mass as it relies on the 
assumption of virial equilibrium ( i]5.3| l. 

The differential mass function for GMC-only 13 CO clumps, 
using the LTE estimate , is shown in Figure [9] Following 
Williams & McKee] h997h . we have fitted the distribution with a 
truncated power law, 



where f3 is the usual power law for the mass distribution 
(dN/dM oc M p ). Guided by the fact that CPROPS can only iden- 
tify peaks which lie well above the extended background emis- 
sion from the GMC, we have fitted the distribution only where 
log(M/Mo) > 2, which we adopt as our completeness limit; this 
corresponds (in the optically thin limit) to a clump with a size twice 
the beam area and an average column density equivalent to 7 + 3<j, 
where 7=40 K km s _1 is the average integrated 13 CO intensity and 
a=0.7 K km s _1 is the RMS noise in the integrated 13 CO map. 
In other words, although we are able to detect much lower mass 
clumps, they would not necessarily be identified as local maxima 
in the cube given our adopted decomposition algorithm. 

The fitted power l aw, weighted by the inverse variance assum- 
ing Poission noise (see lWilliams & McKedfl997l) . has an index of 
[3 — —1.68 ± 0.05 for a bin size of 0.2 dex. Varying the bin size 
between 0.1 and 0.25 dex yields power laws between —1.6 and 
— 1.75, suggesting a slightly larger uncertainty in /3 (0.08). Our de- 
rived value of P lies within the range of values (—1.5 to —1.9) that 
have been previously derived for clumps within molecular clouds 
using the iterative GAUSSCLUMPS decomposition method, which 
models the clumps as 3-dimensional Gaussians in position- velocity 
space (e.g.,|Krameretal. 1998; S imon et al. 2001 ). It is also consis- 



tent with the value of — 1.7±0.3 derived by Mookerie a et alj2004h 
using CLUMPFIND and GAUSSCLUM PS for the RCW 106 cl oud 
as observed in dust continuum. Recently, lReid& Wilson] d2006h rc- 
derived the clump mas s function in RCW 106 using the data from 
Mook eriea et al.l 00041) and found a double power law with a slope 
at the high-mass end of —2.4; we do not see indication for such 
a distribution, but further study using denser gas tracers would be 
valuable. 

The dot-dashed histogram shows the clump mass function pro- 
duced using the ECLUMP method. Since this method partitions all 
of the cloud flux among the peaks identified by CPROPS, the result- 
ing clump sizes and masses are significantly larger. These masses 
should be viewed with great caution, as the clump sizes are deter- 
mined largely by the separation between peaks, which in turn was 
determined largely by our decomposition parameters. Note that in 
Paper I, we also found a relatively flat mass function when apply- 
ing CLUMPFIND to the integrated 13 CO map; this may be due 
in part to the difficulty of assigning extended emission to clumps. 
We believe that the steeper slope found with the CPROPS tech- 
nique should be considered more robust, since extended emission 
has been excluded from the clump assignment process. Nonethe- 
less, the form of the mass function is clearly sensitive to the de- 
compostion method employed, especially for a low-density tracer 
such as 13 CO which shows both compact and extended emission. 



4.3 Size-linewidth and mass-radius correlations 

In Figure [10] we examine, in two sets of four panels each, corre- 
lations among clump properties for the 13 CO and C 18 clumps 
within the GMC. Horizontal and vertical dashed lines represent 
our resolution limits; these are given by Al/=0.34 km s" 1 and 
_R=0.64 pc (the latter based on the effective radius for a 45" Gaus- 
sian). For each correlation we compute the slope of the OLS bi- 
sector ( iFeigelson & Babulll992l) . which treats both variables sym- 
metrically by determining the line which bisects the standard or- 
dinary least-squares fits of X on Y, OLS(X\Y), and Y on X, 
OLS(yjX). The quoted errors represent the range between the 
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Figure 10. (Left 4 panels) Correlation plots of derived clump properties for the GMC-only 13 CO emission. Dotted lines are OLS bisector fits; the power-law 
slope is given in the upper left corner. The range between OLS(X| Y) and OLS(Y \X) is indicated by the quoted errors. Vertical and horizontal dashed lines 
represent the instrumental resolution; the skewed dashed line in the luminosity-radius plot represents the sensitivity limit. (Right 4 panels) Correlation plots of 
derived clump properties for the C 18 emission. 
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Figure 11. (Left 4 panels) Correlation plots of derived clump properties for the GMC-only 13 CO emission, using the ECLUMP method. Dotted and dashed 
lines are as per Fig. 1101 (Right 4 panels) Comparison between original values, using the default decomposition method, with the extrapolated values derived 
by CPROPS. Extrapolation increases all values significantly, but adds substantial scatter. 
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OLS{X\Y) and OLS(Y|X) slopes. The virial and LTE masses 
show a clear correlation (upper left panels) but with considerable 
scatter, with the 13 CO clumps tending to have larger virial masses 
compared to LTE masses and the reverse tending to be true for the 



C O clumps. We discuss these trends further in j 

A weak correlation between linewidth and radius appears to 
exist (upper right panels), especially for the C 18 clumps, al- 
though much of the trend is driven by clumps smaller than the 
nominal resolution limit. The slope is poorly constrained because 
of the large scatter in the data points and the limited range of sizes 
probed. The fitted slope ap pears closer to 1 than to 0.5, the latter 
being the value derived by ISolomon et al. for 273 clouds 

observed in CO. Similar studies of correlations fo r clumps within 
clouds (e.g. jKramer et al]l 996; Simo n et al .1200 lh also show con- 
siderable scatter. Superposition effects can strongly affect mea- 
sured clump linewidths, increa sing them beyond th e actual (3-D) 
linewidth for a given size scale l lOstriker et alfeOOll) . 

A clear correlation between luminosity and radius is found for 
both samples (lower left panels). While such a correlation is unsur- 
prising, the slope of the correlation (7 « 3.2) appears steeper than 
would be expected for a constant average column density, which 
would yield M ~ Mlum oc L oc R . Including the slight non- 
linearity between Mlte and Mlum (Fig.[8} would make the agree- 
ment even poorer. The steepness of the correlation cannot be solely 
due to the sensitivity limit, which follows an L oc R 2 relation 
(dashed line). Excluding clumps based on the resolution and com- 
pleteness limits does not significantly change the observed slope, 
although it could then be marginally consistent with 7=2. The 
constant volume (rather than column) density implied by 7 ~ 3 
could be related to the critical density needed to excite the CO 
lines, which is comparable to the average densities of the clumps 
(riff ~ 10 3 ' 3 cm -3 ). Assuming that gas at much higher or lower 
density is uncommon or poorly traced by CO, one would expect 
most clumps to have densities near the critical density. 

Finally, we note that the luminosity and linewidth show a clear 
correlation (lower right panels), although this relation is not inde- 
pendent of the previous two. 

It is worth considering the sensitivity of these relations to the 
decomposition method employed. The left side of Figurell llshows 
the correlations for the 13 CO emission using the ECLUMP method. 
Here both the virial and LTE masses are considerably larger, pre- 
sumably a result of the larger clump size. There is also an even 
steeper relation of luminosity with radius (7 ~ 4). At first glance 
this result is somewhat surprising, as it implies larger clumps have 
higher densities. However, it appears to reflect the tendency for 
smaller clumps to be found at the edges of the cloud: because of the 
competitive nature of how flux is assigned to peaks in this method, 
clumps have a tendency to 'grow' in size in the inner parts of the 
cloud, with the strongest peaks tending to collect the most mass. 

One obvious concern raised by Fig.QiJJis that roughly half (45 
per cent) the clump radii fall below the nominal resolution limit set 
by the Mopra beam. This is possible because the clumps are iden- 
tified as intensity peaks, and the size of those peaks above some 
threshold may be smaller than a beam width. As noted in i]3.3l 
CPROPS can apply corrections for sensitivity and resolution biases, 
extrapolating the sizes, linewidths, and fluxes to a zero noise level 
and then deconvolving the instrumental resolution. Unfortunately, 
these corrections do not appear to be reliable for the 13 CO data, as 
the extrapolated values are poorly correlated with the raw values, 
and can increase the clump masses by orders of magnitude (Fig- 
ure QT| right). Indeed, the sum of the extrapolated clump masses 
exceeds the total GMC mass by a considerable margin (factor of 
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Figure 12. Contours of peak 13 CO opacity overlaid on 5.8 /im GLIMPSE 
image of the Hll region G333.6— 0.2. Contour levels represent optical 
depths of 3, 6, and 9. 



1.3-1.6). Deconvolving the beam without extrapolation excludes 
many clumps as unresolved, as indicated by the resolution limit in 
Fig. [To] but the resulting correlations between properties are largely 
the same, except that the luminosity-radius correlation has a slope 
closer to the constant column density case (7 ~ 2) due to a de- 
crease in size of the smallest clumps. 

We conclude, in agreement with Sim on etld] 1 2001 ), that the 
classi cal scaling relations for molecular clouds (e.g., Solomon et al.l 
Il987h are not easily revealed by studies of substructure within 
clouds, since CO observations probe a limited range of clump size 
and density. Subsequent studies making use of molecular lines of 
higher critical density or dust extinction, both of which can bet- 
ter isolate high-density peaks, should be able to demonstrate these 
relations more clearly. 



4.4 Regions of highest opacity 

As Fig. [3] indicates, the regions of highest opacity in the map are 
not necessarily the regions that have highest 13 CO integrated in- 
tensities. This is not surprising, as high optical depths generally 
require narrow linewidths as well as significant column density. 
The highest optical depths are found in the northern part of the 
map, in a ring-like structure surrounding the luminous H II region 
G333.6— 0.2. Figure [121 shows a 5.8/im image of this region as ob- 
served by the Spitzer Space Telescope as part of the GLIMPSE 
survey (Benjamin et al. 2003). This H II region has been estimated 
bv lFuiivoshi et al.l d2005l) to be between 1 s and 10 6 yr o ld, based 
on the lack of methanol maser emission ( Caswe 13 119971) and the 
presence of 08V stars. The ring, although clumpy, may represent 
a limb-brightened shell driven by stars within the H II region. At a 
projected distance of ~9' (~9 pc) from G333.6— 0.2, an expansion 
velocity of ~ 10 km s~ 1 would be adequate for the shell to reach its 
current size within 10 6 yr. 

As shown in Fig. [T2] several of the high-opacity clouds ap- 
pear fainter than the background emission in the GLIMPSE im- 
age, confirming that they are cold and dense enough to be seen 
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Figure 13. Dependence of column density and opacity on the adopted abun- 
dance ratio T13/T18. Both curves are drawn assuming T m |->( 13 CO)=5 K 
and r mb (C 18 0)=l K. The discontinuous jump in TV at left is due to the 
assumption of T cx =20 K in the optically thin limit. 

in absorption. Given recent interest in these 'infrared dark clouds' 
as sites of current or future massive star formation dSimon et al.l 
2006; Rathbo rne et al 1 l2006t) . a more detailed study of this com- 
plex would be worthwhile. 



5 DISCUSSION 

5.1 Uncertainties in column density 

One limitation of this analysis is the lack of an independent con- 
straint on the kinetic (and thus excitation) temperature at the res- 
olution of the Mopra data. Pointed 12 CO observations by one of 
us (BM) suggest peak TX( 12 CO) values of up to 20-30 K, which 
translates to 40-60 K (T m b). Such brightnesses are roughly a factor 
of 2 larger than the corresponding peak 12 CO brightness tempera- 
tures seen in the NANTEN data. In the optically thin approxima- 
tion, the column density for a thermalized population increases by 
roughly a factor of ~5 for T cx increasing from 10-100 K. Thus 
we underestimate the column density by at most a factor of ~2 
in the warmer regions. With a new spectrometer that is capable of 
observing 12 CO, 13 CO, and C 18 simultaneously, future Mopra 
observations should be able to better constrain the variation in T ex 
across clouds. Even with matched resolution observations, how- 
ever, there would remain some ambiguity about whether the same 
material was being traced by 12 CO, given its high optical depth. 

An even greater limitation in warmer regions of the cloud is 
the use of the LTE approximation, especially in populating the un- 
observed higher rotational levels. Direct measurements of higher J 
transition s are still required to te st the LTE assumption. In nearby 
globules, Harjunpaa et al. (2004) find that the LTE approximation 
overestimates column densities by 10-30 per cent, due to subther- 
mal excitation of the higher rotational levels. In massive GMCs like 
RCW 106, where overall temperatures are likely to be higher, the 
total column density becomes even more sensitive to the column 
densities in excited rotational states, and thus the uncertainty in ex- 
trapolating from N( J = 0) is likely to become more significant. 
On the other hand, volume densities are also likely to be higher, 
which may help to thermalize the higher J transitions. Future ob- 
servations with submillimetre telescopes such as NANTEN2 will 
be important in determining whether a non-LTE approach is re- 
quired. 

Finally, uncertainties and possible variations in the 



CO/C O abundance ratio will affect the inferred opacity 
and column density values. As noted in ij3.ll about 12 per cent of 
the voxels have 13 CO/C 18 intensity ratios exceeding our adopted 
abundance ratio of 7.4. Unusually high abundances of 13 CO 
relative to other CO is otopomers have been noted in the outskirt s 
of molecular clouds toickmanetai] 1 19791 : lLanger et alj [l980), 
and have been considered to result from chemical fractionation, 
in which isotopic exchange reactions ( 13 C + + 12 CO <-> 12 C + 
+ 13 CO) favo r 13 CO at low temperatures dWatson et ail 1 19761: 
lLangen 1 Il977l) . For the isotopic exchange to be effective, the 
medium must be partially irradiated by UV photons (yielding a 
significant abundance of ionised carbon) whilst still being cold 
enough to suppress the reverse exchange reaction. Both conditions 
could be fulfilled at the edges of a cloud. Figure [13] shows how 
the LTE column density (cf. Eqs. Q] H gjl varies for 13 CO/C ls O 
abundance ratios ranging from 5 to 10, assuming T mb ( 13 CO)=5 K, 
T mb (C 18 0)=l K, and a default T ox =20 K (in the optically thin 
limit). Whereas T13 is very sensitive to the adopted abundance 
ratio, the column density 7V( 13 CO) is much less so, as a result 
of the anti-correlation between r and T ex for a given observed 
brightness temperature. 

We conclude that errors in the column density of a factor of >2 
could quite plausibly result from uncertainties in excitation temper- 
ature, isotopic abundance ratio, and the LTE approximation, with 
the greatest uncertainties in warmer regions or cloud edges, where 
the column density may be underestimated. 

5.2 The PDF of the column density 

The probability distribution function (PDF) of the volume 
density in a turbulent medium is well-described by a log- 
normal distributio n, provided that the medium is approxi- 
mately isothermal Jvk zquez-Semademl 1 19941 ; IPadoan et al.lll997l 
IPassot & Va zquez-Semadeni 1993)- Even if the turbulence is not 
strictly isothermal, the PDF should remain close to log-normal, 
devel oping a quasi-power-law tail jPassot & Vazquez-Semadenil 
1998). However, volume densities are difficult to infer observa- 
tionally, requiring observations of several J transitions to be inter- 
preted within the context of a radiative t ransfer model. In or der to 
facilitate comparison with observati ons, lOstriker et al.1 d200ll) and 
IVazquez-Semadeni & Garcfal d200ll) have investigated the PDF of 
the column density in molecular cloud models. For the isothermal 
case, the column density PDF remains log-normal, at least until the 
line-of-sight depth becomes much larger than the autocorrelatio n 
function of the density field (Vazquez-Semade ni & Garcfa1l200ll) . 
Our observation of a log-normal column density PDF lends sup- 
port to the view that this form of the PDF is a generic consequence 

of interstellar turbulence. 

Recently, lLombardiet"ai]d2006h have investigated the column 
density PDF in the Pipe nebula, based on near-infrared (near-IR) 
extinction towards background stars. They also find an approxi- 
mately log-normal PDF, although there are additional peaks in the 
distribution, some of which may be related to additional clouds 
along the line of sight. Comparing their results with previous 12 CO 
mapping, they find that 12 CO fails to trace very low column den- 
sities (due presumably to UV dissociation) and very high column 
densit ies (due presumably to saturation or depletion). Rid ge et all 
d2006l) also find that near-IR extinction yields a much closer match 
to a log-normal PDF than 12 CO and 13 CO line intensities in Ophi- 
uchus and Perseus. Our ability to recover a log-normal distribution 
from obsevations of CO isotopomers is likely a reflection of three 
effects: (1) we correct for optical depth, whereas previous studies 
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Figure 14. Virial parameter defined in Eq. [8] plotted against LTE mass for (a) 13 CO emission from the GMC; (b) non-GMC 13 CO emission; (c) C 18 
emission. Dashed lines are OLS bisector fits. The number of clumps and slope of each fit are indicated on each plot. Note the scale of the horizontal (mass) 
axis differs among the panels. 



have tended to ignore opacity; (2) our observations are confined 
to largely molecular regions where 13 CO is likely tracing the bulk 
of the gas mass; and (3) the GMC contains many active Hll re- 
gions, where depletion of CO onto grains is probably less impor- 
tant than in nearby dark clouds. Thus, the approach we outline here 
may be useful for studying the column density PDF in clouds near 
the Galactic plane where extinction-based methods face greater dif- 
ficulties. 



5.3 Dynamical state of the clumps 

A commonly use d indicator of the dynami cal state of clumps is the 
virial parameter (Bertoldi & McKee 1992), 

= 5agfl _ Mvir 
a GM ~ M L te ' 

which measures the ratio of kinetic to gravitational energies. For 
a> 1, the clump must be confined by external pressure in order 
to survive as a long-lived entity, whereas a ~ 1 indicates t hat the 
clump is close to self-gravitating. iBertoldi & McKed |j992) found 
that for a > 1 (pressure-confined clumps), the virial parameter 
correlates with cloud mass as approximately a oc M~ 2 ^ 3 . 

Figure[l4]shows plots of a vs. A/lte for the GMC-only 13 CO 
clumps, the non-GMC 13 CO clumps, and the C ls O clumps. All 
show a trend of decreasing a wit h M that is roughly consist ent with 
the slope of -2/3 predicted by IBertoldi & McKed dl992h . Some- 
what unexpectedly, however, the more massive clumps within the 
GMC, particularly the C ls O clumps, show values of a <C 1, and 
thus seem unable to support themselves against gravity. By con- 
trast, the non-GMC clumps, which we identify with foreground or 
background Galactic gas, appear to have lower masses and to be 
either pressure supported or marginally self-gravitating. 

The low values of a, down to ~0.1, are also puzzling given 
that studies of clumps w ithin other clouds t ypically find very 
few clumps with a < 1 jKramer et at]|l996l ; iHever et alj|200ll ; 
Simo n et al. iinni). One possibility is that we are overestimating 
the gas mass by assuming thermalisation of all J levels. However, 
Table Q] suggests that this can at most explain a factor of 2 dis- 
crepancy, and it would be counteracted by our likely underestimate 
of the excitation temperature ( ij5.ll ). A second possibility is that 
we underestimate the virial mass because we measure the size and 



linewidth of each clump only near the emission peak, rather than 
extrapolating to infinite sensitivity. However, such extrapolation in- 
creases Mlum and A/lte as well, and thus appears to have little 
net effect on a. A third possibility is that we overestimate the dis- 
tance, since a oc D~ , Since we adopt the near kinematic distance, 
however, it seems unlikely that the molecular complex can be much 
nearer. While none of these effects seem able to explain the low a 
values alone, a combination of them might be a viable solution. 

If values of a < 1 are not an artifact of our analysis tech- 
nique, then there must be additional terms in the virial equation 
that are counteracting gravity. For instance, there may be a signif- 
icant degree of magnetic support for the clumps, as suggested by 
iBot et alj {2007) for clouds in the Small Magellanic Cloud with 
similarly low values of a. In this scenario, one would still need to 
explain why the a- M relation shows a smooth variation from the 
pressure-confined regime (a > 1) to the magnetically supported 
regime (a < 1), and how a strong degree of magnetic support is 
maintained in the most massive clumps, where ambiopolar diffu- 
sion should be most effective because of lower ionisation fractions 
in more shielded regions of the cloud. In addition, as noted below 
( ij5.4t . the axis ratios of the clumps are not consistent with an oblate 
population, as would be expected in the simplest magnetically sup- 
ported models. 

Alternatively, iDib et alj J2007t) have criticized the use of the 
virial parameter to assess the dynamical state of clumps and cores, 
because of its neglect of the surface terms in the virial equation, 
which represent the exchange of energy between a region and its 
environment. A useful illustration is the notion of turbulent pres- 
sure, which is generally assumed to provide support against grav- 
ity, yet turbulent compression might plausibly act to accelerate col- 
lapse. The net effect of turbulent motions depends sensitively on the 
surface kinetic energy term, which is difficul t to constrain observa- 
tionally. The simulations of IDib et"aL I J2007t) show a oc M 
even for a < 1. Yet many of the simulated clumps with a < 1 are 
not truly bound, i.e. they have positive total energy, and are thus 
likely to be transient. This is especially true of the lower density 
regions that would be probed by CO. In future work, we plan to in- 
vestigate the a-M relation for clumps identified in tracers of higher 
gas densities to test this possibility. 
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Figure 15. Intensity-weighted mean velocity image of the GMC-only 13 CO 
emission, with approximate position angles of the cloud and the velocity 
gradient overlaid. A line of Galactic latitude fe=— 0.5° is also shown. Con- 
tours range from —60 to —40 km s — 1 in increments of 5 km s _1 . 



5.4 GMC velocity gradient and clump elongations 

iBains et all ( 2006) (Paper I) had noted a velocity gradient across the 
cloud of ~0.2 km s" 1 arcmin" 1 by inspection of position-velocity 
slices. We measured the gradient to be 0.13 km s _1 arcmin" 1 by 
fitting a plane to an image of the intensity-weighted first moment 
along the velocity axis. Since l'~l pc at our adopted distance, this 
is consistent with the typical velocity gradients, projected on the 
plane of the s ky, of 0.15 km s~ pc - for mol ecular cloud s in th e 
Galactic ring dKoda et alj|2006h . As noted by iKoda et all d2006l) . 
such gradients are substantially larger than the gradient expected 
from Galactic rotation (0.04 km s _1 pc" 1 ). Indeed, large scale ve- 
locity gradients, inconsistent with Galactic rotation, are com monly 
revea led by principal component analysis of CO data cubes ( Brunt 
2003). These gradients may be an indication of turbulence driven 
on large scales, as microturbulent models have difficulty account- 
ing for such gradients (Brunt 20031). 

The observed velocity gradient (Figure [TBI has a position an- 
gle of — 19° (counterclockwise relative to north) and appears there- 
fore somewhat misaligned with the elongation of the cloud (at 
P.A.~25°); it makes an angle of 64° relative to the Galactic plane, 
whereas the misalignment of the cloud itself relative to the plane is 
only ~20°. Since the cloud has not been mapped completely to its 
edges, the orientation of the velocity gradient is not strongly con- 
strained. Nonetheless, an alignment of the gradient with the Galac- 
tic plane appears unlikely. This is consistent with the analysis of 
IKoda et all HoOrj) , who found that molecular clouds in the inner 
Galaxy are preferentially elongated along the Galactic plane, but 
that their spin axes have random orientations. They argue that this 
elongation cannot be due to rotation and must result from internal 



velocity anisotropy, and thus favor a driving source for turbulence 
related to Galactic rotation. However, it is unclear to what extent 
the tendency for clouds to be elongated along the plane is related to 
internal motions rather than the initial conditions for cloud forma- 
tion in a vertically stratified gas layer. 

Focusing on the clumps within the RCW 106 cloud, with typ- 
ical sizes of ~1 pc and line widths of ~1 km s _1 , the large-scale 
velocity gradient appears to play only a minor role in accounting 
for the observed velocity dispersions. Whether the clumps them- 
selves have significant velocity gradients is unclear; since they are 
only slightly resolved in most cases, reliable fitting is difficult and 
measurement biases would need to be carefully quantified. Of the 
four brightest C 18 clumps, three show gradients of >0.5 km s" 1 
arcmin" 1 when fitting a plane to the first moment image, with po- 
sition angles of —31°, 36°, and 118°. We suspect that clump ve- 
locity gradients will show only a weak correspondence with the 
large-scale velocity gradient, just as the morphological position an- 
gles show little alignment with the cloud as a whole or the Galactic 
plane (Fig. [7}. Thus, clumps and their properties appear to be de- 
fined primarily by turbu lence w i thin th e cloud. 

As pointed out by Rvden (1996), mean apparent axis ratios 
of 0.5-0.6, as fou nd by the present work and previous work (e.g., 
iMvers et al .1 1 1 99 lb are fully consistent with a population of pro- 
late objects, but only marginally consistent with a population of 
oblate objects. Therefore static, magnetically supported models, 
which yield oblate equilibri a unless the magnetic fi eld has a sig- 
nificant toroidal component dFiege & P udritz 2000), are probably 
not applicable for the clumps we observe. Further supp ort for this 
conclusion comes from work by iGammie et alj d2003l) , who ap- 
plied a clump-finding algorithm to 3-D simulations of magneto- 
hydrodynamic turbulence. Their distribution of apparent axis ratios 
peaks around 0.5, and only 10 per cent of their clumps are oblate. 
However, our observations lack the spatial resolution to study the 
elongation of star-forming cores, where the importance of magnetic 
support is most controversial. 



6 CONCLUSIONS 

In this paper we have presented C 18 mapping of the central part 
of the GMC associated with RCW 106, which had been previously 
mapped at similar resolution in 13 CO (Paper I). The 13 CO and 
C ls O data have been used to estimate the optical depth and excita- 
tion temperature as a function of position and velocity, and thereby 
obtain estimates of the column density assuming LTE. In a sepa- 
rate analysis, we have decomposed both 13 CO and C 18 datacubes 
into clumps using an automated decomposition algorithm and esti- 
mated the clump masses using the LTE column density cube. The 
principal results of this study are as follows: 

(i) Correcting for opacity does not have a dramatic effect on the 
overall cloud mass, since regions of high opacity are restricted to 
compact regions. Depending on the adopted excitation temperature, 
assuming negligible opacity can even lead to an overestimate of the 
cloud mass. 

(ii) Correcting the 13 CO emission for opacity does broaden the 
distribution of column densities, which otherwise exhibits a rel- 
atively rapid falloff at the high end. The resulting column density 
PDF is well-approximated as a log-normal function, consistent with 
what is commonly found in numerical simulations of turbulence. 

(iii) The regions of highest opacity are found in the northern part 
of the cloud, in a ring-like structure surrounding the luminous H II 
region G333.6— 0.2. We speculate that this ring is a limb-brightened 
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shell, ~10 pc in radius, driven by the expansion of the H II region. 
The shell material is cold enough to appear in absorption in mid- 
infrared Spitzer images. 

(iv) The distribution of clump masses can be approximated as a 
power law, dN/dM oc M~ 1,7 , consistent with previous studies of 
clumps within molecular clouds. While this result appears sensitive 
to the decomposition method employed, we argue that the adopted 
method is less prone to assigning unrelated extended emission to 
the largest clumps. 

(v) Positive correlations between clump size, linewidth, and 
luminosity are found, albeit with considerable scatter. The 
luminosity-radius correlation appears closer to L oc R 3 than L oc 
R 2 , suggesting most clumps have a volume density close to the 
critical density needed to excite CO emission. The strength and in- 
terpretation of such correlations should be revisited in future work 
with maps tracing higher density gas. 

(vi) The virial parameter a, which measures the ratio of kinetic 
to gravitational energies, is inversely correlated with the LTE clump 
mass. The most massive clumps in the GMC show d<l, indicat- 
ing that another source of support against gravity is required, pos- 
sibly magnetic fields. Alternatively, a may fail to capture the true 
dynamical state of the clumps, due to the neglect of time-varying 
surface terms in the virial theorem. 

(vii) The clumps show mean apparent axis ratios of ~0.5, with 
no preferential alignment along the Galactic plane. On the other 
hand, the GMC as a whole does appear elongated along the plane, 
although it exhibits a large-scale velocity gradient with a different 
orientation. 

Future studies in this series will probe the turbulent proper- 
ties of the cloud using the power spectrum and related techniques, 
and incorporate observations of higher-density tracers, as well as 
the young stellar population, to address directly the relationship 
between molecular cloud properties and star formation. 
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